AD-A162  883 


l# 


National  Research 
Council  Canada 


Conseil  national 
de  recherches  Canada 


m 


SUBSONIC  WALL 
INTERFERENCE  CORRECTIONS 
FOR  HALF-MODEL  TESTS 
USING  SPARSE 
WALL  PRESSURE  JTATA 


[::l 


DTIC 

ELECTE 
JAN  0  7  1986 


by 

M.  Mokry 
National  Aeronautical  Establishment 


OTTAWA 
NOVEMBER  1985 


f  *  •  r  • 


mm 


This  document  has  been  approved 
ioi  public  release  and  sale;  its 
distribution  is  unlimited. 


Canada 


AERONAUTICAL  REPORT 
LR-616 

NRC  NO.  25132 


86 


'd2*ZcJjLj£  * 


*  •  *  •  •  m  *  •  ■  V  *  >  •  j,  ■  at  *  •  *  ■  *  »  •  b  *  |*  *  *  * 


i  %  m 


NATIONAL  AERONAUTICAL  ESTABLISHMENT 
SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 

AERONAUTICAL  REPORTS: 

Aeronautical  Reports  (LR):  Scientific  and  technical  information  pertaining  to  aeronautics  considered 
important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

Mechanical  Engineering  Reports  (MS):  Scientific  and  technical  information  pertaining  to  investigations 
outside  aeronautics  considered  important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

AERONAUTICAL  NOTES  (AN):  Information  less  broad  in  scope  but  nevertheless  of  importance  as  a 
contribution  to  existing  knowledge. 

LABORATORY  TECHNICAL  REPORTS  (LTR):  Information  receiving  limited  distribution  because 
of  preliminary  data,  security  classification,  proprietary,  or  other  reasons. 

Details  on  the  availability  of  these  publications  may  be  obtained  from: 

Publications  Section, 

National  Research  Council  Canada, 

National  Aeronautical  Establishment, 

Bldg.  M-16,  Room  204, 

Montreal  Road, 

Ottawa,  Ontario 
K1A  0R6 

ETABLISSEMENT  AERONAUTIQUE  NATIONAL 
PUBLICATIONS  SCIENTIFIQUES  ET  TECHNIQUES 

RAPPORTS  D ’AERONAUTIQUE 

Rapports  d’aeronautique  (LR):  Informations  scientifiques  et  techniques  touchant  l’aeronautique 
jugees  importantes,  completes  et  durables  en  termes  de  contribution  aux  connaissances  actuelles. 

Rapports  de  genie  mecanique  (MS):  Informations  scientifiques  et  techniques  sur  la  recherche  exteme 
a  l’aeronautique  jugees  importantes,  completes  et  durables  en  termes  de  contribution  aux  connais¬ 
sances  actuelles. 

CAHIERS  D’AERONAUTIQUE  (AN):  Informations  de  moindre  port4e  mais  importantes  en  termes 
d’accroissement  des  connaissances. 

RAPPORTS  TECHNIQUES  DE  LABORATOIRE  (LTR):  Informations  peu  disseminees  pour  des 
raisons  d’usage  secret,  de  droit  de  propriete  ou  autres  ou  parce  qu’elles  constituent  des  donnees 
preliminaires. 

Les  publications  ci-dessus  peuvent  etre  obtenues  a  l’adresse  suivante: 

Section  des  publications 

Conseil  national  de  recherches  Canada 

Etablissement  aeronautique  national 

Im.  M-16,  piece  204 

Chemin  de  Montreal 

Ottawa  (Ontario) 

K1A  0R6 


SUBSONIC  WALL  INTERFERENCE  CORRECTIONS  FOR  HALF-MODEL  TESTS 
USING  SPARSE  WALL  PRESSURE  DATA 


CORRECTIONS  SUBSONIQUES  DES  EFFETS  DE  PAROIS  POUR  DES  ESSAIS 
DE  DEMI-MAQUETTES  A  L’AIDE  DE  DONNEES  DISPERSEES  DE  PRESSION  SUR  LES  PAROIS 


by/par 
M.  Mokry 


Presented  at  the  Euromech  Colloquium  No.  187  Adaptive  Wall  Wind  Tunnels  and  Wall  Interference 
Correction  Methods,  Gottingen,  October  15-17,  1984. 

Presente  au  colioque  Euromech  n°  187  Adaptive  Wall  Wind  Tunnels  and  Wall  Interference  Correction 
Methods,  Gottingen,  15-17  octobre  1984. 


This  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited.  _ 


L.H.  Ohman,  Head/Chef 

High  Speed  Aerodynamics  Laboratory/ 

Laboratoire  d’aerodynamique  a  hautes  vitesses 


G.M.  Lind  berg 
Director/Directeur 


I '  -  55  I  • .  '  '  il*i»  'fr^V-^J CiLi h-tw *wS **  L-ii^  mL -|L<<  -li  t .  +**  » JL  * ■  M, 


I  wfc  Wl  «!^  »,  VUJLI  W*:4L 


SUMMARY 

A  method  is  described  for  correcting  subsonic  wind  tunnel  measurements  on  half¬ 
models  in  ventilated  test  sections,  operated  at  subcritical  flow  conditions  at  the  walls.  For 
perforated  walls,  the  boundary  values  of  the  streamwise  component  of  the  wall  interference 
velocity  are  obtained  from  static  pressures  measured  by  a  few  longitudinal  pressure  tubes 
or  rails  attached  to  the  walls  and  from  the  estimated  farfield  of  the  model  in  free  air.  The 
sparse  boundary  data  is  extended  by  means  of  streamwise  smoothing  and  transverse  inter¬ 
polation.  The  streamwise  velocity  correction  is  derived  from  the  doublet  panel  solution  of 
an  interior  Dirichlet  problem  and  the  transverse  corrections  by  integrating  the  irrotational 
flow  conditions. 

The  evaluated  corrections  to  Mach  number  and  angle  of  attack,  presented  as  contour 
plots  in  the  wing  plane,  provide  insight  into  the  correctability  of  the  test  results.  Examples 
are  given  for  a  transport  aircraft  half-model  tested  in  the  NAE  high  speed  wind  tunnel. 
Applicability  of  the  method  at  supercritical  flow  conditions  at  the  model  is  examined  on 
experimental  and  computational  data  of  a  high  aspect  ratio  wing  in  the  Lockheed  CFWT 
facility. 


RESUME 

On  decrit  une  methode  de  correction  de  mesures  subsoniques  effectuees  en  soufflerie 
sur  des  demi-maquettes  placees  dans  des  veines  d’essai  ventilees,  dans  des  conditions 
d’ecoulement  subcritique  sur  les  parois.  Pour  des  parois  perforees,  les  valeurs  limites  de  la 
composante  longitudinale  de  la  vitesse  de  perturbation  induite  par  les  parois  sont  obtenues 
a  partir  des  pressions  statiques  mesurees  par  quelques  tubes  de  pression  longitudinaux  ou 
par  des  rails  fixes  aux  parois  et  a  partir  du  champ  lointain  estime  du  modele  a  Fair  fi¬ 
bre.  Les  donnees  limites  dispersees  sont  augmentees  au  moyen  d’un  lissage  longitudinal 
et  d’interpolations  transversales.  La  vitesse  longitudinale  est  corrigee  par  resolution  d’un 
probleme  interieur  de  Dirichlet  par  la  methode  des  doublets  et  la  vitesse  transversale  par 
integration  des  conditions  d’ecoulement  irrotationnel. 

Les  corrections  evaluees  du  nombre  de  Mach  et  de  I’angle  d’incidence,  presentees 
sous  forme  de  contours  dans  le  plan  de  l’aile,  renseignent  sur  la  mesure  dans  laquelle  on 
peut  corriger  les  resultats  des  essais.  Des  examples  sont  donnes  pour  une  demi-maquette 
d’avion  de  transport  essayee  dans  la  soufflerie  a  vitesse  elevee  de  l’EAN.  L’applicabilite  de 
la  methode  a  des  conditions  d’ecoulement  supercritiques  sur  la  maquette  est  examinee  a 
partir  de  donnees  expermentales  et  calculees  pour  une  aile  a  allongement  eieve  placee  dans 
la  soufflerie  a  ecoulement  compressible  (CFWT)  de  Lockheed. 
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V  associated  with  volume 
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1.0  INTRODUCTION 
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The  introduction  of  wind  tunnels  with  adaptive  walls,  making  extensive  use  of  mea¬ 
surements  of  flow  parameters  on  control  surfaces,  gave  a  new  impetus  for  the  development 
of  more  rigorous  correction  methods  for  conventional  wind  tunnels.  After  several  decades 
of  testing  and  comparisons  of  theory  with  experimental  results  it  became  evident  that  ap¬ 
plication  of  idealized  boundary  conditions  representing  flow  through  ventilated  walls  was 
inadequate.  Capclier,  Chevallier,  and  Bounio!  gave  a  pioneering  paper  ll]  on  the  subject, 
showing  that  with  the  help  of  measured  wall  pressures  and  the  farfield  representation  of 
the  model  by  singularities,  it  is  possible  to  obtain  subsonic  corrections  as  a  solution  of  the 
interior  Dirichlet  problem,  with  no  reference  to  the  crossflow  properties  of  the  walls. 

For  solid  walls,  a  method  utilizing  wall  pressure  measurements  was  independently 
proposed  by  Ashill  and  Weeks  [2].  The  major  difference  between  the  two  techniques  is  due 
to  the  fact  that  at  a  solid  wall  the  velocity  vector  is  fully  determined  by  static  pressure. 
The  corrections  can  be  obtained  from  Green’s  theorem,  making  the  farfield  representation 
of  the  model  superfluous.  Hackett  et  al.[3]  showed  that  for  solid  walls  the  location  of  model 
singularities  and  their  strengths  can  in  fact  be  deduced  from  wall  pressure  signatures  alone. 
The  advantage  of  not  requiring  the  simulation  of  the  model  farfield  is  of  particular  benefit  in 
cases  where  the  flows  are  complex  and  hence  are  not  easily  modeled,  for  example  those  with 
extensive  separations.  The  method  of  Ref. [2]  is  in  principle  also  applicable  to  ventilated 
tunnel  walls,  however,  a  fast,  reliable  measurement  the  boundary  flow  angles  is  still  an 
unresolved  problem. 

The  “inordinate”  amount  of  experimental  data,  needed  to  specify  (Dirichlet)  boundary 
values  on  bounding  surfaces,  appears  as  an  obstacle  to  the  three-dimensional  interference 
assessment  [4j.  Fortunately,  physics  of  pressure  disturbance  propagation  and  direct  wind 
tunnel  observations  indicate  that  wall  interference  effects  are  in  general  less  severe  in 
three  dimensions  than  those  in  two  dimensions.  Accordingly,  simplifying  assumptions 
can  be  made  to  develop  cost-effective  procedures  that  minimize  the  number  of  boundary 
measurements  needed  to  assess  wall  interference  [5].  Thus  Ref.[l]  suggests  to  use  the 
Fourier  expansion  of  the  sparse  pressure  data  on  the  bounding  surfaces,  showing  that  that 
in  fact  only  a  first  few  terms  are  required  to  achieve  an  acceptable  accuracy.  Sawada  [6] 
devised  a  similar  technique,  giving  a  comprehensive  analysis  of  the  weight  functions  of  the 
Fourier  components.  Mokry  [7]  represented  the  model  by  a  point  disturbance  and  evaluated 
the  corrections  from  the  Fourier  solution  for  the  Dirichlet  problem  in  a  finite-length  circular 
cylinder.  Wall  pressures  measured  by  four  longitudinal  pressure  tubes  were  used  to  define 
the  lowest  order  harmonics,  that  proved  to  be  sufficient  for  defining  the  Mach  number  and 
the  flow  angle  corrections  at  the  model  position.  Conceivably,  this  technique  is  limited  to 
small  models  in  test  sections  permitting  the  placement  of  pressure  tubes  on  a  cylindrical 
surface,  such  as  those  having  circular,  octagonal  or  square  cross-sections.  The  method 
of  Rizk  and  Smithmeyer  [8],  developed  during  the  same  period,  is  applicable  to  general 
rectangular  cross-sections  and  uses  a  more  realistic  representation  of  the  test  model.  Static 
pressures  are  assumed  to  be  known  all  over  the  test  section  boundary,  using  a  combination 
of  sparse  measurements  and  theoretical  estimates.  A  potential  flow  model  of  a  slotted  wall 
test  section,  intended  for  use  with  sparsely  measured  wall  pressues,  has  more  recently  been 


proposed  by  Kemp  [4].  While  this  approach  is  suitable  for  longitudinally  slotted  walls, 
preserving  streamwise  momentum,  it  is  apparently  not  applicable  to  perforated  walls, 
having  viscous  cross-flow  properties. 

The  present  method  is  devised  for  rectangular  tunnel  cross-sections  and  has  many 
common  features  with  that  of  Ref.[8].  For  perforated  test  section  walls,  the  boundary 
pressures  are  assumed  to  be  measured  by  a  number  of  pressure  tubes  or  rails,  Fig.l,  at¬ 
tached  directly  to  the  walls  as  shown  on  schematic  diagrams  in  Fig. 2.  The  rails,  Fig. lb, 
by  the  nature  of  their  design  [9],  are  excellent  static  pressure  devices  for  two-dimensional 
testing,  but  the  simpler  tubes,  Fig. la,  appear  to  be  more  tolerant  to  three-dimensional 
disturbance  fields,  such  as  generated  by  half-models.  An  asymptotic  analysis  [10]  of  tube 
surface  pressures  indicates  that  the  departure  from  the  ambient  static  pressure  (measure¬ 
ment  error),  being  proportional  to  the  product  of  the  tube  rad'us  and  the  streamwise 
derivative  of  normal  velocity,  was  in  our  case  insignificant  and  could  be  ignored.  From  the 
point  of  view  of  applicability  of  the  method,  the  most  important  feature  was  a  relatively 
low  scatter  of  the  measured  pressures,  see  a  typical  example  in  Fig. 5,  allowing  to  interpret 
them  as  “averaged”  pressures  over  the  open  and  closed  portions  of  the  walls.  For  longitudi¬ 
nally  slotted  tunnel  walls,  the  use  of  similar  pressure  devices  is  very  unlikely,  as  the  pressure 
taps  can  directly  be  installed  on  the  slats.  The  interpretation  of  slat  pressures  in  terms 
of  “averaged”  boundary  pressures  is  not  as  simple  [4j,jll],|l2],  but  the  implementation  of 
the  method  is  in  principle  as  for  perforated  walls. 

The  boundary  values  of  the  streamwise  component  of  wall  interference  velocity,  ob¬ 
tained  by  subtracting  the  estimated  pressure  disturbance  of  the  model  in  free  air,  are 
smoothed  and  expanded  over  the  bounding  surfaces  in  transverse  directions  by  means  of 
interpolations.  The  interior  Dirichlet  problem  for  the  streamwise  component  of  wall  inter¬ 
ference  velocity  is  solved  by  the  first-order  doublet  panel  method.  The  transverse  compo¬ 
nents  of  wall  interference  velocity  are  derived  by  integrating  the  irrot duality  conditions 
and  assuming  that  flow  is  parallel  to  the  tunnel  axis  at  a  reference  plane  far  upstream. 


2.0  DESCRIPTION  OF  THE  METHOD 

It  is  stipulated  that  there  exists  a  region  between  the  tunnel  walls  and  the  model 
where  the  flow  can  be  adequately  described  by  the  linearized  potential  equation  * 
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2  d24>  d2<t>  d2<t> 

dX2  4  dY2  +  dZ2 
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(2.1) 


In  accordance  with  the  classical  concept  of  subsonic  wall  interference  [13],  the  disturbance 
velocity  potential  in  the  linearized  flow  region  is  decomposed  into  two  terms: 


<t>  =  <t>F  +  4>w , 


(2.2) 


*  One  of  the  advantages  of  this  method  is  th;  t  it  only  requires  the  linearized  equation 
to  be  valid  near  the  wall  (excluding  the  shear  layers).  This  means  that  the  procedure  may 
be  valid  as  long  as  supersonic  pockets  do  not  extend  to  the  tunnel  walls. 


where  Of  is  the  disturbance  potential  of  the  model  in  free  air  and  <f>w  is  the  disturbance 
potential,  representing  the  effect  of  tunnel  walls.  It  is  further  assumed  that  <t> w  satisfies 
Eq.(2.1)  in  the  entire  tunnel  interior,  including  the  volume  occupied  by  the  model.  The 
derivatives 

d<t>w  d<t>w  d<pw  ln 

ax  ay  dz  K  1 

are  the  interference  velocity  components,  interpreted  as  disturbances  to  unit  stream  veloc¬ 
ity.  Of  particular  interest  are  their  values  at  the  model,  which  determine  the  local  flow 
conditions  to  which  the  tested  model  is  subjected.  Because  of  differentiability  of  harmonic 
functions,  the  interference  velocity  components  inside  the  test  section  again  satisfy 

2  d2u  d2u  d2u 

0  ax2  +  W2  +  ’  (  ’4) 

and  so  on. 

If  the  static  pressure  is  measured  on  the  tunnel  boundary,  the  boundary  values  of  the 
streamwise  component  of  wall  interference  velocity,  u,  can  be  obtained  via  the  linearized 
pressure  coefficient 

/~<  n  In  r  \ 
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From  Eqs.(2.2)  and  (2.3) 

,,  _  1  r  d(t>F  ,  , 

2 Cp  dX ‘  (2'6) 

The  streamwise  derivative  of  4>f  on  the  tunnel  boundary  can  be  estimated  from  measured 
forces  and  model  geometry  as  discussed  in  Appendix  1. 

Provided  that  the  values  of  u  are  known  over  a  closed  test  section  boundary,  Eqs.(2.4) 
and  (2.6)  specify  an  interior  Dirichlet  problem,  whose  solution  is  known  to  exist  and  be 
unique. 

The  transverse  velocity  components,  v  and  iv,  can  similarly  be  determined  from  their 
respective  (Dirichlet)  boundary  values.  If,  however,  only  the  boundary  values  of  u  are 
known,  the  determination  of  v  and  w  via  irrotational  flow  conditions  is  nonunique.  One 
sees  immediately  that  if  t p  =  ip(Y,Z)  is  an  arbitrary  function  satisfying  Laplace’s  equation 

d2<p  d 2\p 

- -f  — —  —  0 

dY 2  dZ 2 

in  the  transverse  Y,  Z  plane,  then  the  sum  <f>  +  <p  again  satisfies  Eqs.(2.1)  and  (2.5). 

In  two  dimensions  (absence  of  coordinate  Y),  the  integration  yields  d<pjdZ  =  constant 
and  hence  the  transverse  component  of  velocity,  w ,  is  determined  from  u  onl.  up  to  a 
constant  [14  .  In  three  dimensions,  the  interference  velocity  components  v  and  w  are 
determined  from  u  to  within  an  arbitrary  solution  of  Laplace’s  equation  in  the  transverse 
plane.  Integrating  the  differentials  dv  and  dw  in  the  streamwise  direction,  we  obtain  [8]: 


v.v  v.v 
,v  v.v.v 
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-(e.v.z)de. 


w(X,Y,Z)  iv{Xr,Y,Z) 


M,r<z)di, 


V  f 


where  Xr  is  the  A'-coordinate  of  the  (upstream)  reference  plane.  Accordingly,  v  and  w 
can  be  determined  from  Eqs.(2.7)  at  any  point  X.Y,  Z  inside  the  test  section,  provided 
that  their  reference-plane  values  are  known.  Flow  angle  measurement  is  a  possibility  (15], 
but  it  is  unlikely  that  it  could  be  done  in  “production”  testing. 

One  of  the  simple  options  of  bypassing  the  measurement  is  to  assume  that  at  a  (dis¬ 
tant)  reference  plane  the  flow  is  parallel  to  tunnel  axis: 

^(X*,Y,Z)  =  0, 

L  (2-8) 

||(Xft,Y,Z)  =  0. 

In  principle,  the  “zero  flow  angle  conditions”  (2.8)  are  justifiable  for  a  distant  reference 
plane  in  a  closed  wall  test  section,  where  the  theoretical  flow  angle  is  known  to  decrease 
with  the  distance  from  the  model  by  an  order  of  magnitude  faster  than  in  free  air.  To 
examine  their  validity  for  a  ventilated  wall  test  section,  we  differentiate  Eqs.(2.8)  with 
respect  to  Y  and  Z  and  substitute  them  in  Eq.(2.1).  This  shows  that 

d2d> 

-£{XR,Y,Z)  =  0, 
or,  in  terms  of  the  linearized  pressure  coefficient  (2.5), 

^(Xfi,Y,Z)  =  0.  (2.9) 

Accordingly,  Eqs.  (2.8)  are  admissible  as  upstream  boundary  conditions  if  the  measured 
streamwise  pressure  curves  display  zero  slopes  at  the  reference  plane. 

Inserting  Eqs. (2.2)  and  (2.8)  in  (2.3),  we  obtain  the  reference-plane  values 


v(xr,y,z)  =  -^J(xft,y,z), 

w(Xr,Y,Z)=  -^(Xr,Y,Z), 


(2.10) 


that  may  be  evaluated  from  the  farfield  of  see  Appendix  1,  and  substituted  in  Eqs. (2. 7). 

Once  the  interior  values  of  components  of  wall  interference  velocity  have  been  com¬ 
puted,  the  wall  interference  corrections  can  be  evaluated  by  standard  procedures  [16].  Thus 
the  Mach  number  correction  is  obtained  as 


M2)Mu , 


(2.11) 


where  M  is  the  reference  Mach  number  (on  which  the  pressure  coefficients  are  based)  and 
Tf(  1-4)  is  the  ratio  of  specific  heats.  The  sideslip  and  angle  of  attack  corrections  (in 
radians)  are  directly  the  values  of  v  and  w  respectively: 
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Usually,  these  corrections  are  evaluated  at  a  representative  model  station  (centre  of  load¬ 
ing),  or  as  averages  over  the  model.  Additionally,  flow  nonuniformities  that  may  be  of 
relevance  for  large  model  dimensions  can  be  expressed  as  spatial  variations  of  Mach  num¬ 
ber  and  flow  angle  corrections  [17]. 


3.0  DOUBLET  PANEL  METHOD 


The  classical  method  for  solving  the  interior  Dirichlet  problem  is  to  construct  the 
solution  in  terms  of  a  double  layer  distribution  on  the  bounding  surface.  To  facilitate  this 
approach,  we  use  the  compressibility  transformation 

x  =  X,  y  =  0Y,  z  =  0Z  (3.1) 

to  reduce  Eq.(2.4)  to  Laplace’s  equation,  V2u  —  0.  Introducing  the  position  vectors 

f0={xo,yo,z0)  and  r=(i,y,2) 

of  the  observation  point  and  the  running  point  respectively,  the  streamwise  component  of 
the  wall  interference  velocity  can  be  expressed  in  terms  of  a  double  layer  distribution 


u[ro)  =  ff  fir) 


d [ 

dn 


1 


4tt]  T0 


dS, 


(3.2) 


where  S  is  the  test  section  boundary,  the  scalar  function  /  of  the  vector  argument  r  is 
the  doublet  density ,  and  d j dn  denotes  the  outward  normal  derivative.  All  differential  and 
integral  operations  are  performed  with  respect  to  the  unsubscripted  coordinates. 

In  the  limit  fa  — ►  r*  €  S,  as  the  observation  point  f0  approaches  a  smooth  surface 
point  Eq.(3.2)  becomes  the  Fredholm  integral  equation  of  the  second  kind  for  the 
doublet  density  18  : 


u(r0  =  -  |/(rfc) 


d 

dn 


1 

4rr  ir* 


--  )  dS. 

T I 


(3.3) 


The  crossed  integral  sign  is  used  to  point  out  that  a  small  area  surrounding  the  point 
r*  is  excluded  from  the  surface  integral.  The  integral  equation  is  solved  numerically, 
by  approximating  the  surface  by  planar  panels  S},  ( j  -  1 ,  ...,JV),  and  using  piecewise 
constant  doublet  densities, 

=  rESj. 

Since 

d  1  _  n  ■  {fk  r) 

dn  |  fk  -  f  | 

is  zero  if  €  Sk ,  Eq.(3.3)  reduces  to 


u{rk) 


1 


fk 


£’>IL 


j  i 
jjk 


i 


dn  \47r|fk  -  rl 


dS. 


(3.4) 


Selecting  as  collocation  points  r*  all  panel  centroids,  we  obtain  a  system  of  linear  algebraic 
equations 


N 


~  Uk, 


k  =  1 . N 


y=i 

for  the  N  unknown  values  of  doublet  densities  f}.  The  matrix  element 


(3.5) 


if  j  =  fc; 


Akj  - 


(3.6) 


is  evaluated  for  a  rectangular  wall  panel  in  closed  form  in  Appendix  2. 

Once  the  doublet  densities  are  known,  the  value  of  the  streamwise  component  of  the 
wall  interference  velocity  at  an  arbitrary  interior  point  f0  is  calculated  from 


(3.7) 


which  is  the  discrete  form  of  integral  (3.2).  The  transverse  components  of  wind  tunnel 
interference  velocity,  v  and  w,  are  obtained  using  Eqs.(2.7)  and  (2.10). 

An  interesting  property  of  the  double  layer  distribution  (3.2)  emerges  when  setting 
/(r)  =  -1  on  the  bounding  surface  S.  Converting  the  surface  integral  to  a  volume  integral 
and  using  the  fact  that  -  l/(47r|r0  -  r|)  is  the  fundamental  solution  for  the  Laplace  operator, 
we  obtain 


fs  dn  \4tt| 

-//MsKbj)" 


ffl 


6(r0  -  r)dV  -l,  r0  G  V, 


where  V  is  the  volume  enclosed  by  the  surface  S  and  6  is  the  Dirac  delta  function.  Con¬ 
versely,  if  u  is  specified  to  be  unity  on  the  smooth  portions  of  bounding  surface  S,  then 
the  corresponding  doublet  density  /  is  minus  one.  This  property  can  be  used  as  a  simple 
numerical  check  for  the  doublet  panel  program. 


4.0  COMPLETION  OF  BOUNDARY  DATA 

Apart  from  the  uncertainty  of  corrections  utilizing  the  classical  interference  concept, 
the  major  source  of  inaccuracy  of  the  method  is  due  to  the  fact  that  the  boundary  values 
at  panel  centroids  have  to  be  constructed  from  the  sparse  data  provided  by  a  few  pressure 
tubes  attached  to  the  walls.  Based  on  some  preliminary  tunnel  flow  calculations  (19], [20], 
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it  was  observed  that  evaluating  the  boundary  values  it  from  Eq.(2.6)  along  the  pressure 
tubes,  smoothing  them,  and  interpolating  in  transverse  directions  gives  a  good  chance  of 
obtaining  a  meaningful  approximation  to  the  entire  boundary  data. 

Four  tubes  located  as  shown  in  Fig. 2a,  allowing  to  capture  spanwise  variations  of 
pressures  in  the  top  and  bottom  planes,  were  used  by  Digney  |20)  and  Khalid  [21]  in 
NAE  half- model  tests.  Using  s>mmetry  about  the  reflection  plane,  du/dy  =  0  at  y  -  0, 
two  tubes  per  wall  (Tl,  T3  and  T2,  T4)  permit  quadratic  interpolation  in  the  spanwise 
direction.  If  the  sidewall  is  solid,  it  acts  as  an  additional  reflection  plane,  and  cubic 
interpolation  can  be  used.  In  the  four-tube  configuration  of  Fig.2a,  the  sidewall  boundary 
values  can  be  approximated  by  interpolating  the  extrapolated  corner  values  at  the  top 
and  bottom  walls.  If  appreciable  blockage  effects  are  expected  (bulky  fuselage),  additional 
sidewall  tubes  may  be  required,  Figs.2b  and  2c.  The  six-tube  configuration  in  Fig.2c  is 
proposed  for  the  next  phase  of  half-model  measurements  in  the  NAE  5  ft  x  5  ft  test  section. 
An  alternative  six-rail  configuration  in  Fig. 2d,  dictated  by  the  elongation  of  the  20  in  x  28 
in  cross  section  of  the  Lockheed-Georgia  Compressible  Flow  Wind  Tunnel  (CFWT),  was 
utilized  by  Pounds  and  Walker  (22}. 

The  completion  of  Dirichlet  boundary  values  at  the  upstream  and  downstream  ends 
of  the  test  section  may  yet  present  another  problem,  since  supplementary  pressure  mea¬ 
surements  are  difficult  to  make  across  the  stream.  Of  course,  if  identical  Cp  values  are 
attained  at  the  endplanes  by  all  tubes,  Cp  may  be  assumed  constant  at  these  planes  and 
the  boundary  values  of  u  evaluated  quite  simply  from  Eq.(2.6).  Otherwise,  interpolations 
of  the  earlier  constructed  edge  values  need  to  be  employed. 


5.0  APPLICATION  OF  THE  METHOD 

Following  a  successful  verification  of  the  method  on  the  simple  test  procedure  decribed 
at  the  end  of  Sect. 3.0,  the  method  was  applied  to  actual  half-model  measurements  in 
the  NAE  High  Speed  Wind  Tunnel.  The  walls  of  the  5  ft  x  5  ft  transonic  test  section 
are  perforated  by  normal  holes  of  0.5  in  diameter,  providing  an  open  area  ratio  of  20%. 
Boundary  pressures  were  measured  by  four  1  in  diameter  pressure  tubes,  Fig. la,  each 
containing  40  static  pressure  orifices.  The  tubes  were  mounted  to  upper  and  lower  walls, 
as  shown  schematically  on  Fig  2a,  with  the  pressure  orifices  facing  the  model  (x,  y  plane). 
The  transverse  coordinates  of  the  pressure  orifices,  displaced  by  the  tube  diameter  from 
the  walls,  are  given  in  Table  1.  The  streamwise  distribution  of  pressure  orifices  will  become 
apparent  in  the  figures  presenting  wall  panelling  and  measured  pressure  distributions. 

The  divisions  of  the  test  section  boundaries  into  72  and  200  rectangular  panels,  ex¬ 
plored  in  the  present  paper,  are  shown  in  Figs. 3  and  4.  This  “open  box”  panelling  excludes 
the  reflection  plane,  to  which  the  model  is  mounted.  The  locations  of  pressure  taps  of  the 
pressure  tubes  T1-T4  are  indicated  by  rows  of  crosses  on  upper  and  lower  boundaries. 
The  measured  static  pressures  along  the  tubes  are  plotted  in  Fig. 5,  for  a  typical  example: 
M  =  0.60  and  a  =  1.951°,  corresponding  to  a  measured  Ci  -  0.689.  The  corrections  AM 
and  Aa  evaluated  in  the  zero-incidence  wing  plane  are  shown  in  Fig.6  (72  panels)  and  Fig. 7 
(200  panels).  The  difference  between  the  contour  plots  of  Figs. 5  and  6  is  not  large,  so  that 


either  of  them  can  be  used  for  a  qualitative  description.  The  Mach  number  correction  is 
observed  to  be  small  and  fairly  uniform  over  the  wing,  which  is  important  from  the  point 
of  view  of  reliability  of  test  data  on  this  most  sensitive  component  of  the  model.  However, 
the  variation  of  AM  over  the  fuselage  indicates  the  presence  of  a  buoyancy  force.  Since  the 
flow  was  accelerated  over  the  model  (a  negative  A  A /  correction  is  required  at  the  nose  and 
a  positive  correction  at  the  tail),  the  measured  drag  force  was  obviously  larger  than  the 
one  that  would  act  on  the  model  in  free  air.  The  angle  of  attack  correction  varies  less  than 
0.1°  over  the  wing  span,  which  is  within  the  limits  of  acceptance.  Downstream  gradients  of 
AM  and  Aa.  corresponding  to  the  drift  of  measured  pressures  towards  negative  values  at 
the  test  section  exit,  see  Fig. 5.  influence  the  model  data  only  indirectly,  via  trailing  shear 
layers.  On  the  whole,  the  wind  tunnel  data  appears  to  be  “correctable”  to  free  air  data. 


Table  1.  NAE  pressure  tubes 


Tube 

y(in) 

z(in) 

Tl 

11.25 

29.00 

T2 

11.25 

-29.00 

T3 

44.25 

29.00 

T4 

44.25 

-29.00 

A  summary  plot  of  AM  and  Aa  versus  the  measured  Ci,  containing  five  different 
incidence  cases,  is  shown  in  Fig. 8.  The  symbols  indicate  correction  values  evaluated  at 
the  centre  of  theoretical  wing  loading.  The  relatively  small  differences  between  the  200 
panel  solutions  (dark  symbols)  and  the  72  panel  solutions  (open  symbols)  suggest  that  the 
coarser  panelling  can  be  used  for  a  quick  assessment  of  corrections  by  an  on-line  computer 
during  wind  tunnel  tests. 

Applicability  of  the  method  at  supercritical  conditions  at  the  model  was  examined 
on  the  experimental  and  computational  data  [23)  and  (24)  of  a  high  aspect  ratio  wing 
“A”,  tested  in  the  Lockheed-Georgia  Compressible  Flow  Wind  Tunnel  (CFWT).  The  ex¬ 
perimental  data  used  here  were  released  to  AGARD  as  test  cases  for  computers  program 
assessment  [25; .  The  walls  of  the  CFW'T  20  in  x  28  in  test  section,  perforated  with  0.25 
in  diameter  holes  slanted  60°  from  the  normal,  were  adjusted  by  moving  outside  shutter 
plates  to  provide  a  4%  open  area  ratio.  Six  pressure  rails,  each  containing  31  pressure 
orifices  and  designed  in  smaller  scale  similar  to  Fig. lb,  were  used  to  measure  static  pres¬ 
sures  at  the  test  section  walls.  The  rails  were  mounted  above  and  below  the  model  as 
shown  schematically  in  Fig. 2d.  The  actual  y  and  z  coordinates  of  the  pressure  orifices 
are  given  in  Table  2.  The  experimental  rail  pressure  distributions,  obtained  at  M  =  0.82 
and  Cl  -  0.53,  are  presented  in  Fig. 9.  Together  with  the  “A"  Wing  geometry  [25],  they 
provide  sufficient  input  data  for  the  evaluation  of  the  corrections  by  the  present  method. 

The  Mach  number  correction  of  Hinson  and  Burdg°s  23  ,  24  .  obtained  by  an  ex¬ 
tended  Bailey-Ballhaus  transonic  flow  code  26  .  is  used  here  for  comparison.  The  method 
is  based  on  matching  the  theoretical  wing  pressures  computed  using  experimental  pres¬ 
sure  boundary  conditions  and  free  air  (unbounded)  conditions.  The  adjustments  of  free 


-9- 


stream  Mach  number  and  of  angle  of  attack,  required  to  establish  agreement  with  the 
pressure-boundary-condition  computation,  define  the  corrections.  For  Wing  “A”,  a  sim¬ 
ple  correction  to  freestream  Mach  number  of  approximately  AM  =  -0.005  was  found  to 
provide  an  acceptable  match  [23], [24]. 


Table  2.  CFWT  pressure  rails 


Tube 

y(in) 

z(in) 

Tl 

6.00 

9.25 

T2 

5.56 

-9.25 

T3 

12.93 

9.25 

T4 

12.88 

-9.25 

T5 

19.68 

9.25 

T6 

19.62 

-9.25 

The  corrections  produced  by  the  present  panel  method  in  Fig.  10  indicate  that  the  strip 
-0.004  >  AM  >  —0.006  extends  over  the  central  portion  of  the  wing,  and  affirms  that 
the  measurement  is  correctable  for  blockage  by  an  adjustment  of  —0.005  in  Mach  number. 
The  angle  of  attack  correction  is  observed  to  be  less  uniform,  however,  no  comparison  is 
available  in  Refs. [23]  or  [24]  to  confirm  our  result  Aft  ~  -0.3°  near  the  centre  of  wing 
load. 


6.0  CONCLUSIONS 

A  linear  method  has  been  developed  for  the  evaluation  wall  interference  corrections 
for  half-model  configurations  at  subcritical  conditions  at  the  ventilated  test  section  walls. 
Input  data  required  for  the  method  are  static  pressures  measured  along  several  longitudinal 
rows  of  pressure  orifices  near  or  on  the  walls,  global  geometrical  parameters  of  the  model, 
and  measured  forces.  Similarly  to  other  methods  based  on  model  representation  and 
measured  pressure  boundary  conditions,  the  crossflow  properties  of  the  walls  do  not  enter 
the  problem  explicitly.  The  technique  is  simple,  straightforward,  and  suitable  for  a  routine 
post-test  assessment  of  corrections  in  test  section  with  perforated  walls. 
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FIG.  1:  CROSS-SECTIONS  OF  NAE  STATIC  PRESSURE  DEVICES 
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FIG.  3b:  TEST  SECTION  PANELLING,  UNFOLDED  (72  PANELS) 
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APPENDIX  Is  FARFIELD  REPRESENTATION  OF  THE  MODEL 


The  calculation  of  boundary  values  of  u  from  Eq.(2.6)  requires,  in  addition  to  the 
measured  pressure  coefficients,  also  knowledge  of  the  streamwise  derivative  of  the  free  air 
disturbance  velocity  potential  4>F- 

The  disturbance  velocity  potential  due  to  a  wing  of  moderate  to  large  aspect  ratio, 
observed  at  a  farfield  point  (X,  Y,  Z),  can  be  represented  by  discrete  singularities  located 
at  points  (A'm,  Ym,  Zm)  distributed  along  the  quarter-chord  line. 

In  accordance  with  the  lifting  line  theory,  the  disturbance  velocity  potential  associated 
with  lift  can  be  approximated  by  a  discrete  distribution  of  horseshoe  vortices 


4>lf{X,Y,Z)  = 


1  V~"  p  Z  -  Zm 

4^  m(Y-Ymy  +  (Z-ZJf 


x 


- 

V  V(x  -  +  WW  -  iy2  + 


The  individual  strengths  Tm  are  obtained  from  the  spanwise  wing  load  distribution.  The 
total  strength  (circulation)  satisfies 


r  =  £r"  =  \AC<- 

m 


where  Cl  is  the  lift  coefficient  of  the  wing  and  A  is  the  reference  area  on  which  it  is  based. 

Similarly,  the  farfield  associated  with  drag  can  be  represented  by  a  discrete  distribution 
of  sources 


<t>?(X,Y,Z)  =  -  ~Tom 

47T  l — * 

m 


_ 1 _ 

VTx  -  xmy  +  p (y  Y~y^(z  -  zm y ' 


The  total  source  strength  satisfies 


m 


and  Co  is  the  (wake)  drag  coefficient  of  the  wing. 

Finally,  the  farfield  associated  with  the  volume  effect  of  the  wing  can  be  represented 
by  a  discrete  distribution  of  doublets 


<t>Vr(X,Y,Z) 


1  y'  v  _ _ x  :  Xm _ _ _ 

4^  m \(x  -  xmy  +  p(y  -  ym)2  +  p{z  -  zm)2p/2  ’ 


v  = 


where 


is  the  total  volume  of  the  wing.  Alternative  expressions  for  the  disturbance  velocity  po¬ 
tential  due  to  wings  can  be  found  in  Ref.}8'. 

The  fuselage  or  nacelle  can  be  modelled  by  sources  and  sinks,  for  example  as  described 
in  Ref .[3] .  The  simplest  case  is  a  Rankine  body,  generated  by  a  source  and  a  sink  of  equal 
strengths  on  a  line  parallel  to  the  stream  direction.  Thus  for  a  source  located  at  (-£,0,0) 
and  a  sink  at  (£,0,0),  the  disturbance  velocity  potential  is  given  by 

4>f[X,Y,Z)  =  —  1  .  ^  _ 1  . 

*  '  ’  A  _  /  /  *  r  tn\n  mir  o  m  T')  ,  /  ir  r*\  o  n  /in 


4tt  [^(X  +  £)2  +  32Y2  -+  32Z2  x/{X  -  £)2  +  32Y2  +  /?2Z2J 

The  half-axes  A  and  B  can  be  determined  by  evaluating  the  stagnation  point  distance  and 
the  maximum  deflection  of  the  dividing  streamline  respectively  |27 . ,  28  .  The  half-axes  A 
and  B  and  the  “eccentricity”  £  are  found  to  be  interrelated  by 

ab2(e  -t  6)  -  (a  e)4 

where 

£  =  £2,  a  =  A2,  and  b  32  B2 . 

Specifying  the  half-axes  to  approximate  the  shape  of  the  body.  £  can  be  evaluated  in  a  few 
iterations  by  the  generalized  Newton’s  method  (29j,  using 


-  0  -  ¥): 


as  the  initial  guess.  The  singularity  strength  n  follows  from 

[A2  -  £2)2 

M  =  T— 


The  farfield  potential,  combining  the  above  effects,  is  then  given  by  the  sum 

4>f  -  <j>F  +  <t>F  +  Q'f  +  ^ F  ■ 


APPENDIX  2:  VELOCITY  INDUCED  BY  A  RECTANGULAR  PANEL 


The  analysis  is  greatly  simplified  by  introducing  local  coordinates  £,77,5-,  with  f  along 
the  (outward)  normal  to  the  panel.  The  streamwise  velocity  component  at  the  point 
s oiWoito’  induced  by  a  rectangular  panel  £i  <  S<  h'fti  <  V  <  <T  ~  0  of  unit  doublet 
density  is 

| - WT{k=WW.-W]'<-'‘didn 


4 *L  L  Kfo  ; ->)j + £  ” 

fo  p _ 1 _ r  _ _ ]<=«■ 

4rr  ivo  -  v)2  +  f2  V(£°  £)2  +  {Vo  -  r?)2  +  J  C=C. 


\tto  -0J  4  (^o  1  ^)2  +  ?2j3/2 

i  [  { 


revealing  familiar  features  of  a  vortex-type  singularity.  The  remaining  integration  is  per¬ 
formed  with  the  help  of  the  indefinite  integral 


dY  1  B  Y 

_ _ _ _ . _ —  _____  |  _  _ 

(V2  *  A2)\/Y*  +  A2  +  B2  AB  [A  v/y 2  +  A2  +  B2 


The  result  is 


uiZo-r/o^o)  =  ~  {  arctan 
47 r  l 


(o-Cl 


Vo  -  V i 


-  arctan 


arctan 


arctan 


\/(£o  -  fl)2  +  [Vo  -  V\)2  +  ?o  ■ 
jo  -  Cl  _ r?o  -  >/2 _ ' 

\/(£o  -  £l)2  +  (Vo  -  Vi)2  + 

go  -  6  _ *7o  -  >?1 _ ' 

\/Uo  -  C2)2  +  {Vo  -  ??l)2  +  fo  ■ 

go  -  6  _ ?7o  -  V2 _ ' 

fo  \/Uo  -  g2)2  +  (ho  ~  V2)2  +  £2 


Considering  the  principal  value  of  the  arctangent  function,  it  is  easily  verified  that 


lim  u  =  <  ^  2  ’ 

fo  — *o±  {  0, 


on  the  panel; 
outside  the  panel, 


in  accordance  with  Eq.(3.6). 

In  Ref.  30 j,  a  similar  formula  is  given  for  the  trapezoidal  doublet  panel,  which  can  be 
used  for  more  general  test  section  geometries. 
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